Early survival in Atlantic salmon is associated with parental genotypes at loci linked to timing of maturation.

Large effects loci often contain genes with critical developmental functions with potentially broad effects across life-stages. However, the life-stage-specific fitness consequences are rarely explored. In Atlantic salmon, variation in two large-effect loci, six6 and vgll3, is linked to age at maturity, and several physiological and behavioural traits in early life. By genotyping the progeny of wild Atlantic salmon that were planted into natural streams with nutrient manipulations, we tested if genetic variation in these loci is associated with survival in early life. We found that higher early life survival was linked to the genotype associated with late maturation in the vgll3, but with early maturation in the six6 locus. These effects were significant in high-nutrient, but not in in low-nutrient streams. The differences in early survival were not explained by additive genetic effects in the offspring generation, but by maternal genotypes in the six6 locus, and by both parents' genotypes in the vgll3 locus. Our results suggest that indirect genetic effects by large-effect loci can be significant determinants of offspring fitness. This study demonstrates an intriguing case of how large-effect loci can exhibit complex fitness associations across life stages in the wild and indicates that predicting evolutionary dynamics is difficult.

Tissue-specific metabolic enzyme levels covary with whole-animal metabolic rates and life-history loci via epistatic effects.

Metabolic rates, including standard (SMR) and maximum (MMR) metabolic rate have often been linked with life-history strategies. Variation in context- and tissue-level metabolism underlying SMR and MMR may thus provide a physiological basis for life-history variation. This raises a hypothesis that tissue-specific metabolism covaries with whole-animal metabolic rates and is genetically linked to life history. In Atlantic salmon (Salmo salar), variation in two loci, vgll3 and six6, affects life history via age-at-maturity as well as MMR. Here, using individuals with known SMR and MMR with different vgll3 and six6 genotype combinations, we measured proxies of mitochondrial density and anaerobic metabolism, i.e. maximal activities of the mitochondrial citrate synthase (CS) and lactate dehydrogenase (LDH) enzymes, in four tissues (heart, intestine, liver, white muscle) across low- and high-food regimes. We found enzymatic activities were related to metabolic rates, mainly SMR, in the intestine and heart. Individual loci were not associated with the enzymatic activities, but we found epistatic effects and genotype-by-environment interactions in CS activity in the heart and epistasis in LDH activity in the intestine. These effects suggest that mitochondrial density and anaerobic capacity in the heart and intestine may partly mediate variation in metabolic rates and life history via age-at-maturity. This article is part of the theme issue 'The evolutionary significance of variation in metabolic rates'.

Standard metabolic rate does not associate with age-at-maturity genotype in juvenile Atlantic salmon.

Atlantic salmon (Salmo salar) is a species with diverse life-history strategies, to which the timing of maturation contributes considerably. Recently, the genome region including the gene vgll3 has gained attention as a locus with a large effect on Atlantic salmon maturation timing, and recent studies on the vgll3 locus in salmon have indicated that its effect might be mediated through body condition and accumulation of adipose tissue. However, the cellular and physiological pathways leading from vgll3 genotype to phenotype are still unknown. Standard metabolic rate is a potentially important trait for resource acquisition and assimilation and we hypothesized that this trait, being a proxy for the maintenance energy expenditure of an individual, could be an important link in the pathway from vgll3 genotype to maturation timing phenotype. As a first step to studying links between vgll3 and the metabolic phenotype of Atlantic salmon, we measured the standard metabolic rate of 150 first-year Atlantic salmon juveniles of both sexes, originating from 14 different families with either late-maturing or early-maturing vgll3 genotypes. No significant difference in mass-adjusted standard metabolic rate was detected between individuals with different vgll3 genotypes, indicating that juvenile salmon of different vgll3 genotypes have similar maintenance energy requirements in the experimental conditions used and that the effects of vgll3 on body condition and maturation are not strongly related to maintenance energy expenditure in either sex at this life stage.

Genetic coupling of life-history and aerobic performance in Atlantic salmon.

A better understanding of the genetic and phenotypic architecture underlying life-history variation is a longstanding aim in biology. Theories suggest energy metabolism determines life-history variation by modulating resource acquisition and allocation trade-offs, but the genetic underpinnings of the relationship and its dependence on ecological conditions have rarely been demonstrated. The strong genetic determination of age-at-maturity by two unlinked genomic regions (vgll3 and six6) makes Atlantic salmon (Salmo salar) an ideal model to address these questions. Using more than 250 juveniles in common garden conditions, we quantified the covariation between metabolic phenotypes-standard and maximum metabolic rates (SMR and MMR), and aerobic scope (AS)-and the life-history genomic regions, and tested if food availability modulates the relationships. We found that the early maturation genotype in vgll3 was associated with higher MMR and consequently AS. Additionally, MMR exhibited physiological epistasis; it was decreased when late maturation genotypes co-occurred in both genomic regions. Contrary to our expectation, the life-history genotypes had no effects on SMR. Furthermore, food availability had no effect on the genetic covariation, suggesting a lack of genotype-by-environment interactions. Our results provide insights on the key organismal processes that link energy use at the juvenile stage to age-at-maturity, indicating potential mechanisms by which metabolism and life-history can coevolve.

Does parental angling selection affect the behavior or metabolism of brown trout parr?

The behavior of organisms can be subject to human-induced selection such as that arising from fishing. Angling is expected to induce mortality on fish with bold and explorative behavior, which are behaviors commonly linked to a high standard metabolic rate. We studied the transgenerational response of brown trout (Salmo trutta) to angling-induced selection by examining the behavior and metabolism of 1-year-old parr between parents that were or were not captured by experimental fly fishing. We performed the angling selection experiment on both a wild and a captive population, and compared the offspring for standard metabolic rate and behavior under predation risk in common garden conditions. Angling had population-specific effects on risk taking and exploration tendency, but no effects on standard metabolic rate. Our study adds to the evidence that angling can induce transgenerational responses on fish personality. However, understanding the mechanisms of divergent responses between the populations requires further study on the selectivity of angling in various conditions.

Resistance is futile: RNA-sequencing reveals differing responses to bat fungal pathogen in Nearctic Myotis lucifugus and Palearctic Myotis myotis.

Resistance and tolerance allow organisms to cope with potentially life-threatening pathogens. Recently introduced pathogens initially induce resistance responses, but natural selection favors the development of tolerance, allowing for a commensal relationship to evolve. Mycosis by Pseudogymnoascus destructans, causing white-nose syndrome (WNS) in Nearctic hibernating bats, has resulted in population declines since 2006. The pathogen, which spread from Europe, has infected species of Palearctic Myotis for a longer period. We compared ecologically relevant responses to the fungal infection in the susceptible Nearctic M. lucifugus and less susceptible Palearctic M. myotis, to uncover factors contributing to survival differences in the two species. Samples were collected from euthermic bats during arousal from hibernation, a naturally occurring phenomenon, during which transcriptional responses are activated. We compared the whole-transcriptome responses in wild bats infected with P. destructans hibernating in their natural habitat. Our results show dramatically different local transcriptional responses to the pathogen between uninfected and infected samples from the two species. Whereas we found 1526 significantly upregulated or downregulated transcripts in infected M. lucifugus, only one transcript was downregulated in M. myotis. The upregulated response pathways in M. lucifugus include immune cell activation and migration, and inflammatory pathways, indicative of an unsuccessful attempt to resist the infection. In contrast, M. myotis appears to tolerate P. destructans infection by not activating a transcriptional response. These host-microbe interactions determine pathology, contributing to WNS susceptibility, or commensalism, promoting tolerance to fungal colonization during hibernation that favors survival.

Association Mapping Based on a Common-Garden Migration Experiment Reveals Candidate Genes for Migration Tendency in Brown Trout.

A better understanding of the environmental and genetic contribution to migratory behavior and the evolution of traits linked to migration is crucial for fish conservation and fisheries management. Up to date, a few genes with unequivocal influence on the adoption of alternative migration strategies have been identified in salmonids. Here, we used a common garden set-up to measure individual migration distances of generally highly polymorphic brown trout Salmo trutta from two populations. Fish from the assumedly resident population showed clearly shorter migration distances than the fish from the assumed migratory population at the ages of 2 and 3 years. By using two alternative analytical pipelines with 22186 and 18264 SNPs obtained through RAD-sequencing, we searched for associations between individual migration distance, and both called genotypes and genotype probabilities. None of the SNPs showed statistically significant individual effects on migration after correction for multiple testing. By choosing a less stringent threshold, defined as an overlap of the top 0.1% SNPs identified by the analytical pipelines, GAPIT and Angsd, we identified eight candidate genes that are potentially linked to individual migration distance. While our results demonstrate large individual and population level differences in migration distances, the detected genetic associations were weak suggesting that migration traits likely have multigenic control.

Comparing RADseq and microsatellites for estimating genetic diversity and relatedness - Implications for brown trout conservation.

The conservation and management of endangered species requires information on their genetic diversity, relatedness and population structure. The main genetic markers applied for these questions are microsatellites and single nucleotide polymorphisms (SNPs), the latter of which remain the more resource demanding approach in most cases. Here, we compare the performance of two approaches, SNPs obtained by restriction-site-associated DNA sequencing (RADseq) and 16 DNA microsatellite loci, for estimating genetic diversity, relatedness and genetic differentiation of three, small, geographically close wild brown trout (Salmo trutta) populations and a regionally used hatchery strain. The genetic differentiation, quantified as F ST, was similar when measured using 16 microsatellites and 4,876 SNPs. Based on both marker types, each brown trout population represented a distinct gene pool with a low level of interbreeding. Analysis of SNPs identified half- and full-siblings with a higher probability than the analysis based on microsatellites, and SNPs outperformed microsatellites in estimating individual-level multilocus heterozygosity. Overall, the results indicated that moderately polymorphic microsatellites and SNPs from RADseq agreed on estimates of population genetic structure in moderately diverged, small populations, but RADseq outperformed microsatellites for applications that required individual-level genotype information, such as quantifying relatedness and individual-level heterozygosity. The results can be applied to other small populations with low or moderate levels of genetic diversity.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis.

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes (e.g., translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr.

Circadian rhythms and environmental disturbances - underexplored interactions.

Biological rhythms control the life of virtually all organisms, impacting numerous aspects ranging from subcellular processes to behaviour. Many studies have shown that changes in abiotic environmental conditions can disturb or entrain circadian (∼24 h) rhythms. These expected changes are so large that they could impose risks to the long-term viability of populations. Climate change is a major global stressor affecting the fitness of animals, partially because it challenges the adaptive associations between endogenous clocks and temperature - consequently, one can posit that a large-scale natural experiment on the plasticity of rhythm-temperature interactions is underway. Further risks are posed by chemical pollution and the depletion of oxygen levels in aquatic environments. Here, we focused our attention on fish, which are at heightened risk of being affected by human influence and are adapted to diverse environments showing predictable changes in light conditions, oxygen saturation and temperature. The examined literature to date suggests an abundance of mechanisms that can lead to interactions between responses to hypoxia, pollutants or pathogens and regulation of endogenous rhythms, but also reveals gaps in our understanding of the plasticity of endogenous rhythms in fish and in how these interactions may be disturbed by human influence and affect natural populations. Here, we summarize research on the molecular mechanisms behind environment-clock interactions as they relate to oxygen variability, temperature and responses to pollutants, and propose ways to address these interactions more conclusively in future studies.

Effect of torpor on host transcriptomic responses to a fungal pathogen in hibernating bats.

Hibernation, the use of prolonged torpor to depress metabolism, is employed by mammals to conserve resources during extended periods of extreme temperatures and/or resource limitation. Mammalian hibernators arouse to euthermy periodically during torpor for reasons that are not well understood, and these arousals may facilitate immune processes. To determine whether arousals enable host responses to pathogens, we used dual RNA-Seq and a paired sampling approach to examine gene expression in a hibernating bat, the little brown myotis (Myotis lucifugus). During torpor, transcript levels differed in only a few genes between uninfected wing tissue and adjacent tissue infected with Pseudogymnoascus destructans, the fungal pathogen that causes white-nose syndrome. Within 70-80 min after emergence from torpor, large changes in gene expression were observed due to local infection, particularly in genes involved in pro-inflammatory host responses to fungal pathogens, but also in many genes involved in immune responses and metabolism. These results support the hypothesis that torpor is a period of relative immune dormancy and arousals allow for local immune responses in infected tissues during hibernation. Host-pathogen interactions were also found to regulate gene expression in the pathogen differently depending on the torpor state of the host. Hibernating species must balance the benefits of energy and water conservation achieved during torpor with the costs of decreased immune competence. Interbout arousals allow hibernators to optimize these, and other, trade-offs during prolonged hibernation by enabling host responses to pathogens within brief, periodic episodes of euthermy.

Cold temperature represses daily rhythms in the liver transcriptome of a stenothermal teleost under decreasing day length.

The climate-change-driven increase in temperature is occurring rapidly and decreasing the predictability of seasonal rhythms at high latitudes. It is therefore urgent to understand how a change in the relationship between photoperiod and temperature can affect ectotherms in these environments. We tested whether temperature affects daily rhythms of transcription in a cold-adapted salmonid using high-throughput RNA sequencing. Arctic char (Salvelinus alpinus) from a subarctic population were reared at a high and a low temperature (15 and 8°C) for 1 month under natural, decreasing day length during late summer. Liver transcriptomes were compared between samples collected in the middle and towards the end of the light period and in the middle of the dark period. Daily variation in transcription was lower in fish from the low temperature compared with strong daily variation in warm-acclimated fish, suggesting that cold temperatures dampen the cycling of transcriptional rhythms under a simultaneously decreasing day length. Different circadian clock genes had divergent expression patterns, responding either by decreased expression or by increased rhythmicity at 15°C compared with 8°C. The results point out mechanisms that can affect the ability of fish to adapt to increasing temperatures caused by climate change.

Pseudogymnoascus destructans transcriptome changes during white-nose syndrome infections.

White nose syndrome (WNS) is caused by the psychrophilic fungus Pseudogymnoascus destructans that can grow in the environment saprotrophically or parasitically by infecting hibernating bats. Infections are pathological in many species of North American bats, disrupting hibernation and causing mortality. To determine what fungal pathways are involved in infection of living tissue, we examined fungal gene expression using RNA-Seq. We compared P. destructans gene expression when grown in culture to that during infection of a North American bat species, Myotis lucifugus, that shows high WNS mortality. Cultured P. destructans was grown at 10 to 14 C and P. destructans growing in vivo was presumably exposed to temperatures ranging from 4 to 8 C during torpor and up to 37 C during periodic arousals. We found that when P. destructans is causing WNS, the most significant differentially expressed genes were involved in heat shock responses, cell wall remodeling, and micronutrient acquisition. These results indicate that this fungal pathogen responds to host-pathogen interactions by regulating gene expression in ways that may contribute to evasion of host responses. Alterations in fungal cell wall structures could allow P. destructans to avoid detection by host pattern recognition receptors and antibody responses. This study has also identified several fungal pathways upregulated during WNS infection that may be candidates for mitigating infection pathology. By identifying host-specific pathogen responses, these observations have important implications for host-pathogen evolutionary relationships in WNS and other fungal diseases.

Molecular Detection of Candidatus Bartonella mayotimonensis in North American Bats.

Candidatus Bartonella mayotimonensis was detected in 2010 from an aortic valve sample of a patient with endocarditis from Iowa, the United States of America. The environmental source of the potentially new endocarditis-causing Bartonella remained elusive. We set out to study the prevalence and diversity of bat-associated Bartonella in North America. During 2015, mist nets and harp traps were used to capture 92 bats belonging to two species: little brown myotis (Myotis lucifugus Le Conte 1831, n = 73) and the gray myotis (M. grisescens A.H. Howell 1909, n = 19) in Kentucky, Michigan, Pennsylvania, and Tennessee. DNA preparations of peripheral blood samples from bats were subjected to a three-marker (gltA, rpoB, and intergenic spacer region [ISR]) multilocus sequence analysis. Sequence-verified gltA-positive PCR amplicons were obtained from nine samples. Three sequences were 99.7-100% identical with the gltA sequence of the Iowa endocarditis patient strain. Analysis of rpoB and ISR sequences demonstrated that one little brown myotis sample from the Upper Peninsula of Michigan contained Bartonella DNA, with 100% sequence identity with the Iowa endocarditis patient strain DNA. It appears possible that bats are a reservoir of Candidatus Bartonella mayotimonensis in North America.

Different Relationship between hsp70 mRNA and hsp70 Levels in the Heat Shock Response of Two Salmonids with Dissimilar Temperature Preference.

The heat shock response (HSR) refers to the rapid production of heat shock proteins (hsps) in response to a sudden increase in temperature. Its regulation by heat shock factors is a good example of how gene expression is transcriptionally regulated by environmental stresses. In contrast, little is known about post-transcriptional regulation of the response. The heat shock response is often used to characterize the temperature tolerance of species with the rationale that whenever the response sets on, a species is approaching its lethal temperature. It has commonly been considered that an increase in hsp mRNA gives an accurate indication that the same happens to the protein level, but this need not be the case. With climate change, understanding the effects of temperature on gene expression of especially polar organisms has become imperative to evaluate how both biodiversity and commercially important species respond, since temperature increases are expected to be largest in polar areas. Here we studied the HSR of two phylogenetically related Arctic species, which differ in their temperature tolerance with Arctic charr having lower maximally tolerated temperature than Atlantic salmon. Arctic charr acclimated to 15°C and exposed to 7°C temperature increase for 30 min showed both an increase in hsp70 mRNA and hsp70 whereas in salmon only hsp70 mRNA increased. Our results indicate that the temperature for transcriptional induction of hsp can be different from the one required for a measurable change in inducible hsp level. The species with lower temperature tolerance, Arctic charr, are experiencing temperature stress already at the higher acclimation temperature, 15°C, as their hsp70 mRNA and hsp70 levels were higher, and they grow less than fish at 8°C (whereas for salmon the opposite is true). Consequently, charr experience more drastic heat shock than salmon. Although further studies are needed to establish the temperature range and length of exposure where hsp mRNA and hsp level are disconnected, the observation suggests that by measuring both hsp mRNA and hsp level, one can evaluate if a species is approaching the higher end of its temperature tolerance, and thus evaluate the vulnerability of an organism to the challenges imposed by elevated water temperature.

The effects of the painkiller diclofenac and hypoxia on gene transcription and antioxidant system in the gills of three-spined stickleback.

Aquatic organisms face multiple stressors in natural ecosystems. More and more often painkillers are detected in surface waters since their prescription has increased worldwide within the last years. Here we examined the effects of the non-steroidal anti-inflammatory drug (NSAID) diclofenac and hypoxia on three-spined sticklebacks (Gasterosteus aculeatus). We exposed sticklebacks to an environmentally relevant concentration of diclofenac (1µg/L) for 14days, to 24h of hypoxia (2.0mg O2/L), and a combination of both. Hypoxia and diclofenac both can be associated with oxidative stress in fish, but it is unclear whether they would act synergistically. Expression analysis of genes related to antioxidant response, hypoxia response, and chemical metabolism in gills showed that diclofenac alone had little effect, while the combination of hypoxia and diclofenac affected transcript levels most, indicating synergistic effects of these stressors. Of the antioxidant enzymes, only superoxide dismutase activity remained unchanged by treatments, while glutathione peroxidase (GPx) was the most affected antioxidant response on both the transcript and activity levels. Our results suggest that diclofenac may lead to suppressed catalase (CAT) activity but increased GPx activity, probably as compensatory mechanism to remove increasing H2O2 in the gills, and that this response is not affected by hypoxia. The activities of lactate dehydrogenase, CAT, and GPx also showed temporal variability during treatments, which can be attributable to tissue-specific circadian rhythms. Our study shows how responses to NSAIDs and hypoxia can interact in fish, suggesting that getting more insight into temporal variation and about the different levels of regulation of environmental responses is necessary in future studies.

Microarray analysis of di-n-butyl phthalate and 17α ethinyl-oestradiol responses in three-spined stickleback testes reveals novel candidate genes for endocrine disruption.

Phthalate esters are plasticizers frequently found in wastewater effluents. Previous studies on phthalates have reported anti-androgenic activity in mammals, causing concerns of their potential effects on the reproduction of aquatic organisms. Another group of environmental endocrine disrupters, steroidal estrogens, are known to inhibit steroid biosynthesis in the gonads, but the effects related to spermatogenesis are not well understood in fish. In this study, three-spined sticklebacks were exposed to di-n-butyl phthalate (DBP) and 17α ethinyl-oestradiol (EE2) at nominal concentrations 35µg/L and 40ng/L, respectively, for four days. The aim of the study was to obtain insight into the acute transcriptional responses putatively associated with endocrine disruption. RNA samples from eight individual male fish per treatment (including controls) were used in microarray analysis, covering the expression of approximately 21,000 genes. In the EE2 treatment the results show transcriptional downregulation of genes associated with steroid biosynthesis pathway and up-regulation of genes involved in pathways related to epidermal growth factor signaling and xenobiotic metabolism. The transcriptional response to DBP was in general weaker than to EE2, but based on enrichment analysis, we suggest adverse effects on retinoid metabolism, creatine kinase activity and cell adhesion. Among the genes showing highest fold changes after DBP treatment compared to control was the teleost fish -specific cytochrome P450 17A2. Overall, this study promotes our understanding on molecular responses to anti-androgens and estrogens in fish testes.

Warm acclimation and oxygen depletion induce species-specific responses in salmonids.

Anthropogenic activities are greatly altering the habitats of animals, whereby fish are already encountering several stressors simultaneously. The purpose of the current study was to investigate the capacity of fish to respond to two different environmental stressors (high temperature and overnight hypoxia) separately and together. We found that acclimation to increased temperature (from 7.7±0.02°C to 14.9±0.05°C) and overnight hypoxia (daily changes from normoxia to 63-67% oxygen saturation), simulating climate change and eutrophication, had both antagonistic and synergistic effects on the capacity of fish to tolerate these stressors. The thermal tolerance of Arctic char (Salvelinus alpinus) and landlocked salmon (Salmo salar m. sebago) increased with warm acclimation by 1.3 and 2.2°C, respectively, but decreased when warm temperature was combined with overnight hypoxia (by 0.2 and 0.4°C, respectively). In contrast, the combination of the stressors more than doubled hypoxia tolerance in salmon and also increased hypoxia tolerance in char by 22%. Salmon had 1.2°C higher thermal tolerance than char, but char tolerated much lower oxygen levels than salmon at a given temperature. The changes in hypoxia tolerance were connected to the responses of the oxygen supply and delivery system. The relative ventricle mass was higher in cold- than in warm-acclimated salmon but the thickness of the compact layer of the ventricle increased with the combination of warm and hypoxia acclimation in both species. Char had also significantly larger hearts and thicker compact layers than salmon. The results illustrate that while fish can have protective responses when encountering a single environmental stressor, the combination of stressors can have unexpected species-specific effects that will influence their survival capacity.

Hypoxia and the pharmaceutical diclofenac influence the circadian responses of three-spined stickleback.

Pollution with low concentrations of pharmaceuticals, especially when combined with low-oxygen conditions (hypoxia), is a threat to aquatic ecosystems worldwide. The non-steroidal anti-inflammatory drug diclofenac is commonly detected in wastewater effluents, and has potential to accumulate in the bile of fish. Diclofenac has been shown to activate aryl hydrocarbon receptor (AHR), which induces transcription in the metabolic enzyme cytochrome P450 1a (cyp1a). Previously, crosstalk has been shown to occur between AHR and hypoxia inducible factor 1 (HIF-1). In addition, both of these transcription factors interact with the proteins regulating circadian (24-h) rhythms in vertebrates. Yet little is known about the significance of these interactions during simultaneous exposure to chemicals and hypoxia in fish in vivo. We exposed wild-caught three-spined sticklebacks (Gasterosteus aculeatus) to diclofenac (1 µg/L, 14 days), hypoxia (2.0 mg/L, up to 24h) and the combination of both. We then analyzed markers of chemical biotransformation (EROD activity, cyp1a and ahr mRNA levels), glycolysis (lactate dehydrogenase (LDH) enzyme activity, ldh and enolase 1a mRNA levels), and the transcription of core circadian clock genes clock and period 1 in liver tissue. Samples were taken at three time points during the light period in order to address disturbances in the circadian variation of metabolic processes. The results show that mRNA levels and LDH activity tended to be lowest before the dark period, but this pattern was disturbed by hypoxia and diclofenac. Diclofenac and hypoxia co-exposure induced EROD activity more strongly than diclofenac exposure alone, while cyp1a mRNA level was increased also by hypoxia and diclofenac alone. LDH activity and mRNA expression showed a clear time-dependent response during hypoxia, which is consistent with the previously suggested decreased accumulation of HIF-1 during the dark period. Furthermore, LDH activity and transcription was disturbed by diclofenac, indicating important effects of environmental pollutants in disturbing natural acclimation. This study demonstrates the need for more studies to understand the potential disturbances in endogenous rhythms caused by environmental pollution in natural populations.

Transcriptional divergence of the duplicated hypoxia-inducible factor alpha genes in zebrafish.

Oxygen availability has been a major force in shaping the physiological evolution of animals. Under reduced oxygen availability (hypoxia) major changes in gene expression are mediated by hypoxia-inducible factors (HIF alphas). Tetrapods have three hif alpha genes, whereas zebrafish (Danio rerio) and other cyprinids have six due to a teleost lineage-specific genome duplication. We studied the transcriptional divergence of the six teleost-specific hif alphas by inspecting the tissue-specific transcription patterns in adult zebrafish and by monitoring the early developmental transcription of normoxia- and hypoxia-grown zebrafish embryos. Overall we observed the highest hif alpha mRNA levels in tissues that are important for hypoxic survival, including the brain, gill and heart. Of the paralogs that have not previously received attention (hif alpha-1A, hif alpha-2B and hif alpha-3B) especially the hif alpha-2B transcription levels suggest functional relevance. The hif alpha-1A/B paralogs that have considerable coding sequence divergence displayed more overall transcriptional divergence than the hif alpha-2A/B paralog pair. The hif alpha-2A/B paralogs that are similarly conserved in coding sequence had a divergent transcription pattern during early development. When zebrafish grown in modest hypoxia were compared to normoxia grown fish, only hif alpha-3A transcription was significantly altered. These results suggest that, in zebrafish, the evolutionary retention of each hif alpha paralog pair has involved unique patterns of coding sequence divergence, adult tissue-specific transcriptional divergence or developmental transcriptional divergence.